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Abstract 



Limits on anomalous WW"f and WW Z couplings are presented from a study 
of WW/WZ — ► evjj events produced in pp collisions at ^fs =1.8 TeV. Results 
from the analysis of data collected using the D0 detector during the 1993- 
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1995 Tevatron collider run at Fermilab are combined with those of an earlier 
study from the 1992-1993 run. A fit to the transverse momentum spectrum of 
the W boson yields direct limits on anomalous WWj and WW Z couplings. 
With the assumption that the WWj and WW Z couplings are equal, we 
obtain -0.34 < A < 0.36 (with An = 0) and -0.43 <Ak< 0.59 (with A = 0) 
at the 95% confidence level for a form-factor scale A = 2.0 TeV. 
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I. INTRODUCTION 



The Tevatron pp collider at Fermilab offers one of the best opportunities to test trilinear 
gauge boson couplings [0-0, which are a direct consequence of the non-Abelian SU{2) x 
U{1) gauge structure of the standard model (SM). The trilinear gauge boson couplings can 
be measured directly from gauge boson pair (diboson) production. Production of WW 
and WZ pairs in pp collisions at ^/s = 1.8 TeV can proceed through s-channel boson 
intermediaries, or a t- or w-channel quark exchange processes as shown in Fig. [I]. There are 
important cancellations between the t- or w-diagrams, which involve only couplings of the 
bosons to fermions, and the s-channel diagrams which contain three-boson couplings. These 
cancellations are essential for making calculations of SM diboson production unitary and 
renormalizable. Since the fermionic couplings of the 7 and W and Z bosons have been well 
tested [|J, we may regard diboson production as primarily a test of the three-boson vertex. 
Production of WW pairs is sensitive to both WW^ and WWZ couplings; WZ production 
is sensitive only to WWZ couplings. 




c) d) 

FIG. 1. Feynman diagrams for WW and WZ production at leading order, a) and c): t- and 
u-channel quark exchange diagrams; b) and d): s-channel diagrams with three-boson couplings. 

A generalized effective Lagrangian has been developed to describe the couplings of three 
gauge bosons |J. The Lorentz- invariant effective Lagrangian for the gauge boson self- 
interactions contains fourteen dimensionless coupling parameters, Ay, Ky, g\ Ay, Kv, g\, 
and g\ {V = Z or 7), seven for WWZ interactions and another seven for WW'y interac- 
tions, and two overall couplings, gww-y = — e and gwwz = — e cot 9w, where e and 6w are 
the positron charge and the weak mixing angle. The couplings Ay and Ky conserve charge 
C and parity P. The couplings g\ are odd under CP and C, g\ are odd under C and P, 
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and ky and Ay are odd under CP and P. To first order in the SM (tree level), all of the 
couplings vanish except g\ and Ky {g\ — gf — /c 7 — Kz — 1)- For real photons, gauge 
invariance in electromagnetic interactions does not allow deviations of gj, gj, and gl from 
their SM values of 1, 0, and 0, respectively. The CP- violating WW'j couplings A 7 and k 7 
are tightly constrained by measurements of the neutron electric dipole moment |J. In the 
present study, we assume C, P and CP symmetries are conserved, reducing the independent 
coupling parameters to k 7 , k z , A 7 , \z and gf . 

Cross sections for gauge boson pair production increase for couplings with non-SM values, 
because the cancellation between the t- and w-channel diagrams and the s-channel diagrams 
is destroyed. This can yield large cross sections at high energies, eventually violating tree- 
level unitarity. A consistent description therefore requires anomalous couplings with a form 
factor that causes them to vanish at very high energies. We will use dipole form factors, 
e.g., Ay(s) = Ay/ (1 + s/A 2 ) 2 , where s is the square of the invariant mass of the gauge-boson 
pair. Given a form-factor scale A, the anomalous-coupling parameters are restricted by S- 
matrix unitarity. Assuming that the independent coupling parameters are k = k 7 = Kz and 
A = A 7 = \ z , tree-level unitarity is satisfied if A < [6.88/((k - l) 2 + 2A 2 ] 1/4 TeV gfj. The 
experimental limits on anomalous couplings can be compared with the bounds derived from 
^-matrix unitarity, and constrain the trilinear gauge-boson couplings only if the limits are 
more stringent than the bounds from unitarity for any given value of A. 

For both WW and WZ production processes, the effect of anomalous values of Ay on 
the helicity amplitudes is enhanced for large s. On the other hand, terms containing Any 
(= k v — 1) grow as V§ in the WZ production process, and as s in the WW production 
process. Limits on Any from the study of WW production are therefore expected to be 
tighter than those from WZ production. 

Since anomalous couplings contribute only via s-channel photon or W or Z boson in- 
termediaries, their effects are expected mainly in the region of small vector boson rapidi- 
ties, and the transverse momentum distribution of the vector boson is therefore partic- 
ularly sensitive to anomalous trilinear gauge-boson couplings. This is demonstrated in 
Fig. |2], which shows the distribution of the W boson transverse momentum p?f in simu- 
lated pp — > WW + X — > eujj + X events for anomalous trilinear gauge boson couplings, 
using a dipole form factor with a scale A = 1.5 TeV, and with the couplings for WW'j and 
WW Z assumed to be equal. 

Trilinear gauge-boson couplings can therefore be measured by comparing the shapes of 
the px distributions of the final state gauge bosons with theoretical predictions. Even if 
the background is much larger than the expected gauge-boson pair production signal as is 
the case for the WW/WZ — > eujj process, limits on anomalous couplings can still be set 
using a kinematic region where the effects of anomalous trilinear gauge boson couplings are 
expected to dominate. 

Trilinear gauge-boson couplings have been studied in several experiments. WWj cou- 
plings have been studied in pp collisions by the UA2 ||, CDF ||, and D0 |To| , [Tll collabo- 



rations using W"f events. The UA2 results are based on data taken during the 1988-1990 
CERN pp collider run at y^i = 630 GeV with an integrated luminosity of 13 pb _1 and the 
CDF and D0 data are from the 1992-1993 and 1993-1995 Fermilab pp runs at = 1.8 
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FIG. 2. The pj( spectrum of generated pp — * WW —>■ eujj events with SM couplings and two 
examples of anomalous couplings. 



TeV. WW Z couplings together with the WW'y couplings have been studied by the CDF 



and D0 collaborations using W boson pair production in the dilepton decay modes [[TT| [T3 
and WW/WZ production in the single-lepton modes jnUnj-ffH] . Experiments at the CERN 
LEP Collider have recently reported results of similar studies ||17|| . 

In this report, we present a detailed description of previously summarized work ]TB| on 
WW and WZ production with one W boson decaying into an electron (or a positron) and 
an antineutrino (or a neutrino) and a second W or Z boson decaying into two jets ||TU . 
Due to the limitation in jet-energy resolution, the hadronic decay of aW boson can not be 
differentiated from that of a Z boson. This analysis is based on the data collected during 
the 1993-1995 Tevatron collider run at Fermilab. From the observed candidate events and 
background estimates, 95% confidence level (C.L.) limits are set on the anomalous trilinear 
gauge boson couplings. The results are combined with those from the 1992-1993 data to 
provide the final limits on the couplings from the D0 analysis. 

Brief summaries of the detector and the multilevel trigger and data acquisition systems 
are presented in Sections |I| and |II[ Sections [TV], [V| and |VI] describe our particle identification 



methods, the data sample, and event selection criteria. Sections |V11| and |V11J| are devoted 
to detection efficiency and background estimates. Results and conclusions are presented in 
Sections [X and 1X1. 



II. THE D0 DETECTOR 



The D0 detector [pOjl , illustrated in Fig. || is a general-purpose detector designed for the 
study of proton-antiproton collisions at ^Js = 1.8 TeV and is located at the D0 interaction 
region of the Tevatron ring at Fermilab. 

The innermost part of the detector consists of a set of tracking chambers that surround 
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Muon Chambers 




FIG. 3. Cutaway view of the D0 detector 

the beam pipe. There is no central magnetic field and jets are measured using a compact 
set of calorimeters positioned outside the tracking volume. To identify muons, an additional 
set of tracking chambers is located outside the calorimeter, with a measurement of muon 
momentum provided through magnetized iron toroids placed between the first two muon- 
tracking layers. 

The full detector is about 13 m high x 11 m wide x 17 m long, with a total weight of 
about 5500 tons. The Tevatron beam pipe passes through the center of the detector, while 
the Main Ring beam pipe passes through the upper portion of the calorimetry, approximately 
2 m above the Tevatron beam pipe. The coordinate system used in D0 is right-handed, 
with the z-axis pointing along the direction of the proton beam (southward) and the y-axis 
pointing up. The polar angle 9 = is along the proton beam direction, and the azimuthal 
angle = along the eastward direction. Instead of 9, we often use the pseudorapidity, rj = 
— ln[tan(#/2)]. This quantity approximates the true rapidity y = 1/2 hi[(E + p z )/(E — p z )], 
when the rest mass is much smaller than the total energy. 

A. Central Detector 

The tracking chambers and a transition radiation detector make up the central detector 
(CD). The main purpose of the CD is to measure the trajectories of charged particles and 
determine the z position of the interaction vertex. This information can be used to determine 
whether an electromagnetic energy cluster in the calorimeter is caused by an electron or by 
a photon. Additional information such as the number of tracks and the ionization energy 
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along the track (dE/dx) can be used to determine whether a track is caused by one or several 
closely spaced charged particles, such as a photon conversion. 

The CD consists of four separate subsystems: the vertex drift chamber (VTX), the 
transition radiation detector (TRD), the central drift chamber (CDC), and two forward 
drift chambers (FDC). The full set of CD detectors fits within the inner cylindrical aperture 
of the calorimeters in a volume of radius r = 78 cm and length / = 270 cm. The system 
provides charged-particle tracking over the region \r]\ < 3.2. The trajectories of charged 
particles are measured with a resolution of 2.5 mrad in <fi and 28 mrad in 9. From these 
measurements, the position of the interaction vertex along the z direction is determined 
with a resolution of 6 mm. 

The VTX is the innermost tracking chamber in the D0 detector, occupying the region 
r = 3.7 cm to 16.2 cm. It is made of three mechanically independent concentric layers of 
cells parallel to the beam pipe. The innermost layer has sixteen cells while the outer two 
layers have thirty-two cells each. 

The TRD occupies the space between the VTX and the CDC; it extends from r = 17.5 
cm to 49 cm. The TRD consists of three separate units, each containing a radiator (393 foils 
of 18 /im thick polypropylene in a volume filled with nitrogen gas) and an X-ray detection 
chamber filled with Xe gas. The TRD information is not used in this analysis. 

The CDC is a cylindrical drift chamber, 184 cm along z, located between r = 49.5 and 
r = 74.5 cm, and provides coverage for \q\ < 1.2. It is made up of four concentric rings of 32 
azimuthal cells per ring. Each cell contains seven sense wires (staggered by 200 jum relative 
to each other to help resolve left-right ambiguities), and two delay lines. The r<p position 
of a hit is determined via the drift time measured for the hit wire and the z position of a 
hit is measured using inductive delay lines embedded in the module-walls of the sense wire 
planes. 

The FDC consists of two sets of drift chambers located at the ends of the CDC. They 
perform the same function as the CDC, but for 1.4 < \rj\ < 3.1. Each FDC package consists 
of three separate chambers: a $ module, whose sense wires are radial and measure the 
<p coordinate, sandwiched between a pair of 6 modules whose sense wires measure the 9 
coordinate. 

B. Calorimeters 

The D0 calorimeters are sampling calorimeters, with liquid argon as the sensitive ioniza- 
tion medium. The primary absorber material is depleted uranium, with copper and stainless 
steel used in the outer regions. There are three separate units, each contained in separate 
cryostats: the Central Calorimeter (CC), the North End Calorimeter (ECN), and the South 
End Calorimeter (ECS). The readout cells are arranged in a pseudo-projective geometry 
pointing to the interaction region. 

The calorimeters are subdivided in depth into three distinct types of modules: electro- 
magnetic sections (EM) with relatively thin uranium absorber plates, fine-hadronic sections 
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(FH) with thicker uranium plates, and coarse-hadronic sections (CH) with thick copper or 
stainless steel plates. There are four separate layers for the EM modules in both the CC 
and EC that are readout separately. The first two layers are 2 radiation lengths thick in 
the CC and 0.3 and 2.6 radiation lengths thick in the EC, and measure the initial longi- 
tudinal shower development, where photons and 7T°s differ somewhat on a statistical basis. 
The third layer spans the region of maximum EM shower energy deposition and the fourth 
completes the EM coverage of approximately 20 total radiation lengths. The fine-hadronic 
modules are typically segmented into three or four layers. Typical transverse sizes of towers 
in both EM and hadronic modules are Arj = 0.1 and A0 = 27r/64 ps 0.1. The third section 
of the EM modules is segmented twice as finely in both 77 and to provide more precise 
determination of centroids of EM showers. 

The CC has a length of 2.6 m, covering the pseudorapidity region |?7| < 1.2, and consists 
of three concentric cylindrical rings. There are 32 EM modules in the inner ring, 16 FH 
modules in the surrounding ring, and 16 CH modules in the outer ring. The EM, FH and 
CH module boundaries are rotated with respect to each other so as to prevent having more 
than one intermodular gap intercepting a trajectory from the origin of the detector. 

The two end calorimeters (ECN and ECS) are mirror-images, and contain four types of 
modules. To avoid the dead spaces in a multi- module design, there is just a single large EM 
module and one inner hadronic (IH) module. Outside the EM and IH, there are concentric 
rings of 16 middle and outer hadronic modules (MH and OH). The azimuthal boundaries 
of the MH and OH modules are also offset to prevent cracks through which particles could 
penetrate the calorimeter. This makes the DO detector almost completely hermetic and pro- 
vides an accurate measurement of missing transverse energy. Due to increase in background 
and loss of tracking efficiency for \rj\ > 2.5, electron and photon candidates are restricted to 
1.5 < \rj\ < 2.5 in the EC. 

In the transition region between the CC and EC (0.8 < \r]\ < 1.4), there is a large 
amount of uninstrumented material in the form of cryostat walls, stiffening rings, and module 
endplates. To correct for energy deposited in the uninstrumented material, we use two 
segmented (0.1 x 0.1 in rj x 0) arrays of scintillation counters, called intercryostat detectors. 
In addition, separate single-cell structures called "massless gaps" are mounted on the end 
plates of the CC-FH modules and on the front plates of EC-MH and EC-OH modules, and 
are used to correct showers in this region of the detector. 

The Main Ring beam pipe passes through the outer layers of the CC, ECN and ECS. 
Beam losses from the Main-Ring cause energy deposition in the calorimeter that can bias 
the energy measurement. The data acquisition system either stops recording data during 
periods of Main-Ring activity near the DO detector, or flags such events. 

C. Muon Detectors 

The DO muon detector is designed to identify muons and to determine their trajectories 
and momenta. It is located outside of the calorimeter, and is divided in two subsystems: 
the Wide Angle Muon Spectrometer and the Small Angle Muon Spectrometer. Since the 
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calorimeter is thick enough to absorb most of the debris from electromagnetic and hadronic 
showers, muons can be identified with great confidence. The muon system is not used in 
this analysis, and is therefore not discussed any further. 

III. MULTILEVEL TRIGGER AND DATA ACQUISITION SYSTEMS 

The D0 trigger system is a multilayer hierarchical system. Increasingly complex tests 
are applied to the data at each successive stage to reduce background. 

The first stage, called Level (LO), consists of two scintillator arrays mounted on the front 
surfaces of the EC cryostats, perpendicular to the beam direction. Each array covers a partial 
region of pseudorapidity for 1.9 < \r)\ < 4.3, with nearly complete coverage over the range 
2.2 < \r)\ < 3.9. The LO system is used to detect the occurrence of an inelastic pp collision, 
and serves as the luminosity monitor for the experiment. In addition, it provides fast 
information on the ^-coordinate of the primary collision vertex, by measuring the difference 
in arrival time between particles hitting the north and south LO arrays; this is used in making 
preliminary trigger decisions. A slower, more accurate measurement of the position of the 
interaction vertex, and an indication of the possible occurrence of multiple interactions, are 
also made available for subsequent trigger decisions. The LO trigger is ~ 99% efficient for 
non-diffractive inelastic collisions. The output rate from LO is on the order of 150 kHz at a 
typical luminosity of 1.6 x 10 31 cm -2 s -1 . 

The next stage of the trigger is called Level 1 (LI). It combines the results from individual 
LI components into a set of global decisions that command the readout of the digitization 
crates. It also interacts with the Level 2 trigger (L2). Most of the LI components, such as 
the calorimeter triggers and the muon triggers, operate within the 3.5 jis interval between 
beam crossings, so that all events are examined. However, other components, such as the 
TRD trigger and several components of the calorimeter and muon triggers, called Level 1.5 
trigger (LI. 5), can require more time. The goal of the LI trigger is to reduce the event rate 
to 100-200 Hz. The primary input for the LI trigger consists of 256 trigger terms, each of 
which corresponds to a single bit, indicating that some specific requirement is met. These 
256 terms are reduced to a set of 32 LI trigger bits by a two-dimensional AND-OR logic 
network. An event is said to pass LI if at least one of these 32 bits is set. The LI trigger also 
uses information based on Main Ring activity. To prevent saturation of the trigger system 
by processes with large cross sections, such as QCD multijet production, any particular 
contributor to the LI trigger can be prescaled. 

The LI calorimeter trigger covers the region up to \rj\ < 4.0 in trigger towers of 0.2 x 0.2 
in rj — <j) space. These towers are subdivided longitudinally into electromagnetic and hadronic 
trigger sectors. The output of the LI calorimeter trigger corresponds to the transverse energy 
deposited in these sectors and towers. 

For the 1993-1995 collider run, an LI. 5 trigger for the calorimeter was implemented 
using the LI calorimeter trigger data and filters based on neighbor sums and ratios of the 
EM and total transverse energies. 

When an event satisfies the LI trigger, the data are passed on the D0 data acquisition 
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pathways to a farm of 48 parallel microprocessors, which serve as event builders as well as 
the L2 trigger system. The L2 system collects the digitized data from all elements of the 
detector and trigger blocks for events that successfully pass Level 1. It applies sophisticated 
algorithms to the data to reduce the event rate to about 2 Hz before passing the accepted 
events on to the host computer for monitoring and recording. The data for a specific event 
are sent over parallel paths to memory modules in specific selected nodes. The accepted 
data are collected and formatted in final form in the nodes, and the L2 filter algorithms are 
then executed. 

The L2 filtering process in each node is built around a series of filter tools. Each tool 
has a specific function related to the identification of a type of particle or event character- 
istic. There are tools to recognize jets, muons, calorimeter EM clusters, tracks associated 
with calorimeter clusters, (sum of transverse energies of jets), and ffir (imbalance in 

transverse energy). Other tools recognize specific noise or background conditions. There are 
128 L2 filters available. If all of the L2 requirements (for at least one of these 128 filters) 
are satisfied, the event is said to pass L2 and it is temporarily stored on disk before being 
transferred to an 8 mm magnetic tape. 

Once an event is passed by an L2 node, it is transmitted to the host cluster, where it is 
received by the data logger, a program running on one of the host computers. This program 
and others associated with it are responsible for receiving data from the L2 system and 
copying it to magnetic tape, while performing all necessary bookkeeping tasks (e.g., time 
stamping, recording the run number, an event number, etc.). Part of the data is sent to an 
event pool for online monitoring. 

A. Electron Trigger 

To trigger on electrons, LI requires the transverse energy in the EM section of a trigger 
tower to be above a programmable threshold. The L2 electron algorithm then uses the full 
segmentation of the EM calorimeter to identify electron showers. Using the trigger towers 
that are above threshold at LI as seeds, the algorithm forms clusters that include all cells in 
the four EM layers and the first FH layer in a region of Arj x A0 = 0.3 x 0.3, centered on 
the tower with the highest Et- The longitudinal and transverse energy profile of the cluster 
must satisfy the following requirements: (i) the fraction of the cluster energy in the EM 
section (the EM fraction) must be above a threshold, which depends on energy and detector 
position; and (ii) the difference between the energy depositions in two regions of the third 
EM layer, covering A77 x A0 = 0.25 x 0.25 and 0.15 x 0.15, and centered on the cell with 
the highest Et, must be within a window that depends on the total cluster energy. 

B. Jet Trigger 

The LI jet triggers require the sum of the transverse energy in the EM and FH sections 
of a trigger tower (A77 x A0 = 0.2 x 0.2) to be above a programmable threshold. The L2 
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jet algorithm begins with an i?<r-ordered list of towers that are above threshold at LI. At 
L2, a jet is formed by placing a cone of given radius TZ, where TZ = \/ Arj 2 + A0 2 , around 
the seed tower from LI. If another seed tower lies within the jet cone, it is passed over and 
not allowed to seed a new jet. The summed E T in all of the towers included in the jet cone 
defines the jet E T . If any two jets overlap, then the towers in the overlap region are added 
into the jet candidate that is formed first. To filter out events, requirements on quantities 
such as the minimum transverse energy of a jet, the minimum transverse size of a jet, the 
minimum number of jets, and the pseudorapidity of jets, can be imposed at this point. 

C. Missing Transverse Energy Trigger 

Rare and interesting physics processes often involve production of weakly interacting 
particles such as neutrinos. These particles usually can not be detected directly. However, 
assuming momentum conservation in a collision allows the momenta of such particles to be 
inferred from the vector sum of the momenta of the observed particles. Since the energy 
flow near the beamline is largely undetected, such calculations are realistic only in the plane 
transverse to the beam. The negative of the vector sum of the momenta of the detected 
particles is referred to as missing E T and denoted by I^t'i it is used as an indicator of the 
presence of weakly interacting particles. At L2, If T is computed using the vector sum of 
all calorimeter and intercryostat detector cell energies with respect to the z position of the 
interaction vertex, which is determined from the timing of the hits in the LO counters. 

IV. PARTICLE IDENTIFICATION 
A. Electron 

Electrons and photons are identified by the properties of the shower in the calorimeter. 
The algorithm loops over all EM towers (At] x A0 = 0.1 x 0.1) with energy E > 50 MeV, 
and connects the neighboring tower with the next highest energy. The cluster energy is then 
defined as the sum of the energies of the EM towers and the energies in the corresponding first 
FH layer. The ratio of the energy in the EM cluster to the total energy (EM energy summed 
with the corresponding hadronic layers), defined as the EM fraction, is used to discriminate 
electrons and photons from hadronic showers. A cluster must pass the following criteria to 
be an electron/photon candidate: (i) the EM fraction must be greater than 90% and (ii) 
at least 40% of the energy must be contained in a single 0.1 x 0.1 tower. To distinguish 
electrons from photons, we search for a track in the central detector that extrapolates to 
the EM cluster from the primary interaction vertex within a window of |A^| < 0.1, and 
|A0| < 0.1. If one or more tracks are found, the object is classified as an electron candidate. 
Otherwise, it is classified as a photon candidate. 
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1. Selection Requirements 

The spatial development of EM showers is quite different from that of hadronic showers 
and the shower shape information can be used to differentiate electrons and photons from 
hadrons. The following variables are used for final electron selection: 

(i) Electromagnetic energy fraction. This quantity is based on the observation that electrons 
deposit almost all of their energy in the EM section of the calorimeter, while hadron jets 
are far more penetrating (typically only 10% of their energy is deposited in the EM section 
of the calorimeter). It is defined as the ratio of EM energy to the total shower energy. 
Electrons are required to have at least 95% of their total energy in the EM calorimeter. 
This requirement loses only about 1% of all electrons. 

(ii) Covariance matrix (if-matrix) \ 2 ■ The shape of any shower can be characterized by the 
fraction of the cluster energy deposited in each layer and tower of the calorimeter. These 
fractions are correlated, i.e., an electron shower deposits energies according to the expected 
transverse and longitudinal shapes of an EM shower and a hadron shower following the 
typical development of a hadronic shower. To obtain good discrimination against hadrons, 
we use a covariance matrix technique. The observables in this method are the fractional 
energies in Layers 1,2, and 4 of the EM sector and the fractional energy in each cell of a 6 x 6 
array of cells in Layer 3 centered on the most energetic tower in the EM cluster. To take 
account of the dependence of the shower shape on energy and on the position of the primary 
interaction vertex, we use the logarithm of the shower energy and the z-position of the event 
vertex as the remaining input observables. The event vertex is determined by extrapolating 
CDC tracks to the z axis, and for more than one possibility, the vertex associated with the 
highest number of tracks is chosen as the event vertex. Using these 41 variables, covariance 
matrices are constructed for each of the 37 detector towers (at different values of 77) based 
on Monte Carlo generated electrons. The Monte Carlo showers are tuned to make them 
agree with our test beam measurements of the shower shapes. The 41 observables for any 
given shower can be compared with the parameters of the appropriate covariance matrix to 
define a x 2 , which is to be be less than 100 for electron candidates in the CC and less than 
200 for the EC. This requirement loses about 5% of all true electrons. 

(iii) Isolation. The decay electron from a W boson should not be close to any other object in 
the event. This is quantified by the isolation fraction. If £'(0.4) is the energy deposited in all 
calorimeter cells within the cone 1Z < 0.4 around the direction of the electron, and EM(0.2) 
is the energy deposited in only the EM calorimeter in the cone 1Z < 0.2, the isolation variable 
is then defined as the ratio J = [£(0.4) - EM(0.2)]/EM(0.2). The requirement J < 0.1 
loses only 3% of the electrons from W boson decays. 

(iv) Track-match significance. An important source of background for electrons is the photon 
from the decay of 7r° or rj mesons. Such photons do not produce tracks, but their trajecto- 
ries can overlap with those of nearby charged particles, thereby simulating electrons. This 
background can be reduced by demanding a good spatial match between the energy cluster 
in the calorimeter and nearby charged tracks. The significance S of the mismatch between 
these quantities is given by S = [(A0/<5a^) 2 + (Az/Saz) 2 ] 1 ^ 2 , where A0 is the azimuthal 
mismatch, Az the mismatch along the beam axis, and the 5 are the resolutions of these vari- 
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ables. This form for S is appropriate for the central calorimeter. For the end calorimeter, r 
replaces z. Requiring S < 5 accepts 95(78)% of the CC(EC) electrons reconstructed in the 
central tracker. 

(v) Track-in-road. All electrons from W — > ev decays are required to have a partially 
reconstructed track along the trajectory between the energy cluster in the calorimeter and 
the interaction vertex. This requirement is found to reject 16(14)% of CC(EC) electrons 
from W boson decay. 

In our analysis, we combine the above quantities to form the electron identification 
criteria. A summary of the selection requirements and their acceptance efficiencies is listed 
in Table |( See Sect|VTI]). 

TABLE I. Electron selection requirements and their acceptance efficiencies for W — > ev events. 



Selection 




cc 




EC 


requirement 




e 




e 


.ff-matrix x 2 


< 100 


0.946±0.005 


< 200 


0.950±0.008 


EM fraction 


> 0.95 


0.991±0.003 


> 0.95 


0.987±0.006 


Isolation 


< 0.10 


0.970±0.004 


< 0.10 


0.976±0.007 


Track match 


< 5 


0.948±0.005 


< 5 


0.776±0.012 


Track-in-road 




0.835±0.009 




0.858±0.006 



2. Electromagnetic Energy Corrections 

The energy scales of the calorimeters were originally set through calibration in a test- 
beam. However, due to differences in conditions between the test beam and the D0 envi- 
ronment, additional corrections had to be implemented. 

The EM energy scales for the calorimeters were determined by comparing the measured 
masses of 7r° — > 77, J/ip — > ee, and Z — > ee to their known values. If the electron energy 
measured in the calorimeter and the true energy are related by E meas = aE true + 5, the 
measured and true mass values are, to first order, related by m meas = am tTue + 8f, where 
the calculable variable / reflects the topology of the decay. To determine a and 5, we fit the 



Monte Carlo prediction to the observed resonances, with a and S as free parameters pi 

h0.03 

-0.21 



The values of a and 5 are found to be a = 0.9533 ± 0.0008 and 5 = -0.16tnS? GeV for the 



CC and a = 0.952 ± 0.002 and 5 = -0.1 ± 0.7 GeV for the EC. 



3. Energy Resolution 



The relative energy resolution for electrons and photons in the CC is expressed by the 
empirical relation = C 2 +J^+^, where E and Et are the energy and transverse energy 
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of the incident electron/photon, C is a constant term from uncertainties in calibration, S 
reflects the sampling fluctuation of the liquid argon calorimeter, and iV corresponds to a 
contribution from noise. For the EC, the in the relation is replaced by E. The sampling 
and noise terms are based on results from the test beam. The noise term measured at the 
test beam agrees with the one obtained in the collider environment (based on the width of 
pedestal distributions). The constant term is tuned to match the mass resolution of both 
observed and simulated Z — > ee events. Table |TJ lists these parameters. 

TABLE II. Parameters for describing the energy resolution of electrons and photons. 



Quantity 


CC 


EC 


C 


0.017 


0.009 


S (VGeV) 


0.14 


0.157 


N (GeV) 


0.49 


1.140 



B. Jets 

In our analysis, jets are reconstructed using a fixed-cone algorithm with radius TZ = 
yj Arj' 2 + A<f)' 2 = 0.5. The algorithm forms preclusters of contiguous cells using a radius of 
7^-preciuster = 0.3 centered on the tower with highest E T . Only towers with E T > 1 GeV are 
included in preclusters. These preclusters serve as the starting points for jet reconstruction. 
An E'r-weighted center of gravity is then formed using the E? of all towers within a radius 
TZ of the center of the cluster, and the process is repeated until the jet becomes stable. A 
jet must have E T > 8 GeV. If two jets share energy, they are combined or split, based on 
the fraction of overlapping energy relative to the Et of the lesser jet. If this shared fraction 
exceeds 50%, the jets are combined. 

Although the TZ = 0.3 cone algorithm is more efficient for jet finding than our larger 
cone size, which leads to undesired merging of jets for high-p^ W or Z bosons, the relatively 
large uncertainties in the measurement of jet-energy for the TZ = 0.3 cones negate their 
advantage, and we therefore choose to use the TZ = 0.5 cone algorithm for our studies. 

1. Selection Requirements 

To remove jets produced by cosmic rays, calorimeter noise, and interactions in the Main 
Ring, we developed a set of requirements based on Monte Carlo studies of jets in such 
environments and on data on noise taken with and without colliding beams. The variables 
used are: 

(i) Electromagnetic energy fraction (em/). As for electrons, this quantity is defined as the 
fraction of the total energy deposited in the electromagnetic section of the calorimeter. A 
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requirement on this quantity removes electrons, photons and false jets from the jet sample. 
Electrons and photons typically have a high EM fraction. False jets are caused mainly 
by background from the Main Ring or by noisy or "hot" cells, and therefore generally 
do not contain energy in the EM section, thereby yielding very low EM fractions. Jets 
with 0.05 < emf < 0.95 are defined as acceptable in this analysis. The efficiency of this 
requirement is 99.9% at Et = 20 GeV and decreases to 99.6% at 100 GeV. 

(ii) Hot cell energy fraction (hcf). The hcf is defined as the ratio of the energy in the cell 
of second highest Et to that of the cell with highest Et within a jet. A requirement on 
this quantity is imposed to remove events with a large amount of noise in the calorimeter. 
Hot cells can appear when a discharge occurs between electrodes within a cell; often this 
does not affect neighboring cells. In this case, hcf is small, which signals a problem, since 
the hcf for a jet should not be small because the energy is expected to be distributed over 
cells. If most of the energy is concentrated in a single cell, it is very likely to be a false jet 
reconstructed from discharge noise. For good jets, hcf is found to be greater than 0.1. The 
efficiency of this requirement is 97.3% at E T = 20 GeV and decreases to 96.9% at 100 GeV. 

(iii) Coarse hadronic energy fraction (chf). This quantity is defined as the fraction of jet 
energy deposited in the coarse hadronic section of the calorimeter. The Main Ring at D0 
passes through the CH modules, and any energy deposition related to the Main Ring will 
be concentrated in this section of the calorimeter. Such jets tend to have more than 40% of 
their energy in the CH region, while standard jets have less than 10% of their energy in this 
section of the calorimeter. All acceptable jets are therefore required to have chf < 0.4. The 
efficiency of this requirement is 99.6% at E T = 20 GeV and decreases to 99.3% at 100 GeV. 

2. Hadronic Energy Corrections 

Since the measured jet energy is usually not equal to the energy of the original parton 
that formed the jet, corrections are needed to minimize any systematic bias. Jet energy 
response affected by non-uniformities in the calorimeter, non-linearities in the response to 
hadrons, emission of particles outside of the 71 = 0.5 cone (often referred to as out-of-cone 
showering), noise due to the radioactivity of uranium, and energy overlap from the products 
of soft interactions of spectator partons within the proton and the antiproton ( "underlying 
event"). The first two effects are estimated using a method called Missing- Et Projection 
Fraction (MPF) p2f. 

The MPF method is based on events that contain a single isolated EM cluster (due 
to a photon or a jet that fragmented mostly into neutral mesons), and one hadronic jet 
located opposite in 0, and no other objects in the event. It is assumed that such events 
do not have energetic neutrinos so that any missing transverse energy can be attributed 
to a mismeasurement of the hadronic jet. The EM-cluster energy is corrected using the 
electromagnetic energy corrections described above. Projecting the corrected ]$t along the 
jet axis determines corrections to the jet energy. This correction is averaged over many 
events in the sample to obtain a correction as a function of jet Et, T], and electromagnetic 
content of the jet. The hadronic energy correction is 20% at E = 20 GeV and 15% at 
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E = 100 GeV, and gradually approaches 10% at high E. 

The impact of out-of-cone showering is estimated using Monte Carlo jet events. Effects 
due to the underlying events and uranium noise are determined in separate studies using 
minimum-bias event data. (Minimum-bias data corresponds to inclusive inelastic collisions 
collected using only the L0 trigger.) 



3. Energy Resolution 



The jet energy resolution has been studied by examining momentum balance in dijet 



events [2jJ. The formula used for parametrizing the relative jet energy resolution is ( 



C 2 + 



+ Table [III] shows the values of the parameters for different r\ regions of the 
calorimeter. 



TABLE III. Jet energy resolution for different regions of the calorimeter. 



rj Region C S (VGcV) N (GeV) 

\rj\ < 0.5 0.00±0.01 0.81±0.02 7.07±0.09 

0.5 < \rj\ < 1.0 0.00±0.01 0.91±0.02 6.92±0.09 

1.0 < M < 1.5 0.05±0.01 1.45±0.02 0.00±1.40 

1.5 < \r]\ < 2.0 0.00±0.01 0.48±0.07 8.15±0.21 

2.0 < hi < 3.0 0.01±0.58 1.64±0.13 3.15±2.50 



C. Neutrinos: Missing Transverse Energy 



The presence of neutrinos in an event is inferred from the Ifa- In this analysis we assume 
that the 1$t in each candidate event corresponds to the neutrino from the decay W — > ev. 



1. Missing 



The missing transverse energy in the calorimeter is defined as I^t = (-^T^+^Ty) 1 ^ 2 ? where 
$t x = - Ei Ei sin(^) cos(0 i ) - Ej and E Ty = - Et E { sin(^) sin(^) - Ej &E 3 y . The 
first sum (over i) is over all cells in the calorimeters, intercryostat detectors and massless 
gaps (see Sec. [11 B| ) . The second sum (over j) is over the E T corrections applied to all 
electrons and jets in the event. This can be used to estimate the transverse momentum of 
any neutrinos in an event that does not contain muons, which deposit only a small portion 
of their energy in the calorimeter. The total missing Et is missing Et from the calorimeter 
corrected for the transverse momenta of any observed muon tracks. Since this analysis does 
not use muons, we will refer to the I$t based on the calorimeters as the true Ifc. 
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2. Resolution in I$t 



For an ideal calorimeter, the magnitude of the components of the Et vector would 
sum to zero for events with no true source of I$t- However, detector noise and energy 
resolution in the measurement of jets, photons, and electrons contribute to the I$t- In 
addition, non- uniform response in the detector also results in Ifo. The Tfir resolution for our 
candidate events is parameterized as a = 1.08GeV + 0.019(2 Et), and is based on studies of 
minimum-bias data ||23|| . The X) Et used in the parameterization is quite reasonable because 
the greater the total amount of transverse energy in the event, the larger the possibility for 
its mismeasurement. 



V. DATA SAMPLE 



The analysis of the WW/WZ — ► evjj process is based on data taken during the 1993- 
1995 Tevatron Collider run (called Run lb). The L0 trigger is used to check the presence of 
an inelastic collision, but is not included in the trigger conditions for V^-boson data. This 
was done to allow studies of diffractive iy-boson production. Our analysis uses the collected 
W —>■ ev data sample, with the L0 trigger requirement imposed offline. The LI trigger used 
in this analysis (called the EM1_1_HIGH trigger) requires the presence of an electromagnetic 
trigger tower with Et > 10 GeV. The LI. 5 trigger then requires the LI trigger tower to have 
Et > 15 GeV and checks that the electromagnetic fraction is greater than 85%. The L2 
component of the trigger (called the EM1_EISTRKCC_MS trigger) requires an isolated electron 
candidate with Et > 20 GeV that has a shower shape consistent with that of an electron 
and E T >15 GeV. 

Additional conditions are imposed on the data to further reduce background. Triggers 
that occur at the times when a proton bunch in the Main Ring passes through the detector 
are not used in this analysis. Similarly, triggers that occur during the first 0.4 seconds of the 
2.4-second antiproton production cycle are rejected. Data taken during periods when the 
data acquisition system or the detector sub-systems malfunctioned are also discarded. With 
these trigger requirements, the integrated luminosity of the data sample is estimated to be 
82.3 ± 4.4 pb _1 [24]. The efficiency and turn-on of the L2 trigger are described in Ref. [25 . 
The trigger efficiency for signal is (98.1 ± 1.9)%. 

Data samples that satisfy two other L2 triggers, the EM1_ELE_M0N and ELE_1_M0N triggers, 
are used for background studies. These triggers select events that have an electron candidate 
with Et > 20 GeV and E T > 16 GeV, respectively. The electron candidates in these samples 
must pass the standard shower-shape requirements, but not the isolation requirement. These 
triggers use the same LI and LI. 5 conditions as the trigger used for signal. 
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VI. EVENT SELECTION 



WW/WZ — > eujj candidates are selected by searching for events with an isolated high- 
er electron, large 1$t, and at least two high- Et jets. Electrons in the candidate sample 
must be in the \rj\ < 1.1 region but away from the boundaries between calorimeter modules 
in (A0 > 0.01), or within the region 1.5 < |?7| < 2.5. Jets in the candidate sample must 
be in the region \r]\ < 2.5. 

The W — > ev decay is defined through the presence of only a single isolated electron 
with E T > 25 GeV and I$t > 25 GeV in the event. The transverse mass of the electron and 
neutrino (Ifir) system is required to be Mt > 40 GeV/c 2 , where Mt = {2E t J^t[^ — cos(0 e — 
j- 1 / 2 . The requirement on the electron Et is sufficiently high to provide an efficiency 
that is independent of Et (the hardware threshold of 20 GeV). Requiring only one electron 
reduces background from Z — > ee production. The requirements on Jf T and M T reduce the 
background contribution from misidentified electrons. 

The W/Z — > jj decay is defined by requiring at least two jets with E T > 20 GeV and an 
invariant mass of the two-jet system consistent with that of the W or Z boson (50 < Mjj < 
110 GeV/c 2 ). The dijet invariant mass (Mjj) is calculated via Mjj = {2E T 1 E T 2 [cosh(i]ji — 
i]j 2 ) — cos(0ji — 0j2)]} 1//2 - If there are more than two jets in the event, the two jets with the 
highest dijet invariant mass are chosen to represent the W (or Z) decay. 

The difference between the pt values of the ev and the two-jet systems is used to reduce 
backgrounds. For WW or WZ production, the Pt(gv) — PrQj) distribution should be 
peaked near zero and have a symmetric Gaussian shape, with the width of the Gaussian 
distribution determined primarily by the jet energy resolution. On the other hand, for 
background such as tt production (see Sec. |V111| ), the distribution should be broader and 
asymmetric (shifted to positive values) due to additional 6-quark jets in the events. Our 
analysis therefore requires \pT(ev) —pT(jj)\ < 40 GeV/c. 

The data satisfying the above selection criteria yield 399 events. Figure ^a shows a scatter 
plot of pt{cv) vs pr{jj) for candidate events that satisfy the two-jet mass requirement. The 
width of the band reflects both the resolution and the true spread in the pt values. Figure [|b 
shows a scatter plot of pr(e^) vs Mjj without the imposition of the two-jet mass requirement. 

VII. DETECTION EFFICIENCY 

A. Electron Selection Efficiency 

The efficiency of electron selection is studied using the Z — > ee event sample from the 
1993-1995 Tevatron collider run using the EM2_EIS_HI trigger. Z — > ee events were selected 
at LI and LI. 5 by requiring two EM towers with Et > 7 GeV at LI, and at least one tower 
with Et > 12 GeV with more than 85% of its energy in the EM section of the calorimeter. At 
L2, the trigger required two electron candidates with E T > 20 GeV that satisfied electron 
shower-shape and isolation requirements. To select an unbiased sample of electrons, we 
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use events in which one of the electrons passes the tag quality requirements: EM fraction 
> 0.90, Isolation < 0.15, ^/-matrix x 2 < 100(200) for CC(EC), and track-match significance 
< 10. The second electron in the event is then assumed to be unbiased. If both electrons 
pass the tag requirements, the event contributes twice to the sample. The efficiency of a 
selection requirement for electrons is given by e = (e s — £&/&)/(l — /&), where e s is the 
efficiency measured in the signal region, is the efficiency measured in the background 
region, and /& is the ratio of the number of background events in the signal region to the 
total number of events in the signal region. The signal region is defined as the region of 
the Z mass peak (86 < m ee < 96 GeV/c 2 ), and the background regions are defined as 
61 < m ee < 71 GeV/c 2 and 111 < m ee < 121 GeV/c 2 . We determine /& in the region 
of the signal using an average of the number of events in the background regions. The 
systematic uncertainties on the efficiencies are estimated from a comparison with efficiencies 
obtained using an alternative method that fits the invariant mass spectrum of two electrons 
to the sum of a Breit-Wigner form convoluted with a Gaussian and a linear dependence for 
the background. Efficiencies from the two methods agree within their uncertainties. The 
track-in-road efficiency is estimated in a similar manner, except that EM clusters with no 
matching track are included as unbiased electrons in the sample. Table | summarizes the 
electron efficiencies. Although these efficiencies are based mainly on Z events with few jets, 
the corrections for > 2 jets are small. 
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B. W/Z — > jj Selection Efficiency 

The ty/Z — > jj selection efficiency is estimated using Monte Carlo WW/WZ —>■ evjj 
events generated with the isajet and pythia |27[ programs, followed by a detailed 



simulation of the D0 detector, and parametrized as a function of pY^ . Figure |5] shows 
the W/Z — > jj detection efficiency e(W — > jj) calculated as the ratio of events after the 
imposition of the two-jet selection requirements relative to the initial number of events. At 
low p T , the detection efficiency is artificially elevated due to the presence of additional jets 
from initial- and final-state gluon radiation (ISR/FSR) that are mislabeled as being decays 
of W or Z bosons. The decrease in the efficiency at high is due to the merging of the 
two jets from a W or a Z boson. The results obtained from isajet are used to estimate the 
efficiencies for identifying the WW/WZ process. 
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FIG. 5. Efficiency for W — ► jj selection as a function of . The decrease in the efficiency at 
high pt is due to the merging of the two jets from the decay of a IF boson. 

The estimated W/Z — > jj efficiency is affected by the jet energy scale, the accuracy 
of the ISR/FSR simulation, the accuracy of the parton fragmentation mechanism, and the 
statistics of the Monte Carlo samples. 

The energy-scale correction has an uncertainty that decreases from 5% at jet Ej- = 20 
GeV to 2% at 80 GeV, and then increases to 5% at 350 GeV. The effect of this uncertainty 
has been studied by recalculating the efficiency with the jet energy scale changed by one 
standard deviation. The largest relative change in the accepted number of events is found 
to be 3%. 

To estimate the uncertainty due to the accuracy of the ISR/FSR simulation and of the 
parton fragmentation mechanism, we use the W/Z — > jj efficiency based on Monte Carlo 
samples generated with pythia. The efficiency obtained using isajet is lower than that for 
pythia, but by less than 10%. We use the efficiencies from isajet because they provide 
smaller yields of WW/WZ events and therefore weaker limits on anomalous couplings. We 
define one-half of the largest difference in isajet/pythia efficiency estimations (5%) as the 
systematic uncertainty attributable to the choice of event generator. 



C. Overall Selection Efficiency 

The overall detection efficiency for WW/WZ — > evjj events assuming SM couplings is 
calculated using two MC methods, coupled with electron-selection and trigger efficiencies 



22 



measured from data. The first MC method uses the isajet event generator followed by a 
detailed simulation of the D0 detector. The second MC method uses the event generator of 
Ref. [0 and a fast simulation program to characterize the response of the detector, isajet 
used the CTEQ2L |28| parton distribution functions to simulate 2500 WW — *> evjj events 
and 1000 WZ — > evjj events with SM couplings. The event selection efficiency for for 
the WW — > evjj signal is estimated as eww = (13.4 ± 0.8)%, and ewz = (15.7 ± 1.4)% 
for the WZ — > evjj signal, where the errors are statistical. The combined efficiency for 
WW/WZ — ► evjj is given by [e ww ■ a ■ B(WW — ► evjj) + e wz ■ o ■ B(WZ — > evjj)] /[a ■ 
B(WW — > evjj) + a ■ B(WZ — ► evjj)] = (13.7 ±0.7)%, where the theoretical cross sections 
of 9.5 pb for WW and 2.5 pb for WZ production |29|, and the W and Z boson branching 
fractions from the Particle Data Group f|, are used in the calculation (a-B(WW — > evjj) = 
1.38 ± 0.05 pb and a ■ B{WZ -> evjj) = 0.188 ± 0.006 pb). 

For the fast simulation, we generated over 30,000 events, with approximately four times 
more for WW production than WZ production, reflecting the sizes of their expected produc- 
tion cross sections. The overall detection efficiencies for the SM couplings were calculated 
as [14.7 ± 0.2(stat) ± 1.2(sys)]% for WW -> evjj and [14.6 ± 0.4(stat) ± l.l(sys)]% for 
WZ — > evjj. The 7.8% systematic uncertainty includes statistics of the fast MC (1%), 
efficiency of trigger and electron identification (1%), Ifir smearing and modeling of the px of 
the WW/WZ system (5%), difference in W — > jj detection efficiencies from the two event 
generators (5%), and the effect of the jet energy scale (3%). The combined efficiency is 
[14.7 ± 0.2(stat) ± 1.2(sys)]%. The combined efficiency estimated using the fast simulation 
is consistent with the value obtained using isajet. 

D. Expected Number of Signal Events 

Using the fast detector simulation and the cross section times branching ratio from the 
event generator of Ref. § (a ■ B(WW -> evjj) = 1.26 ± 0.18 pb, and a ■ B{WZ -> evjj) 
= 0.18 ± 0.03 pb), we estimate the number of expected WW/WZ — > evjj events to be 
17.5 ± 3.0 (15.3 ± 3.0 WW events and 2.2 ± 0.5 WZ events), with the uncertainty (17.1%) 
given by the sum in quadrature of the uncertainty in the efficiency, the uncertainty in the 
luminosity (5.4%), and that in the NLO calculation (14%). 

VIII. BACKGROUND 

The sources of background to the WW/WZ — > evjj process can be divided into two 
categories. The first is instrumental background due to misidentified or mismeasured parti- 
cles, and the other is inherent irreducible background consisting of physical processes with 
the same signature as the events of interest. 
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A. Instrumental Background 



The major source of instrumental background is QCD multijet production in which one 
of the jets showers (mainly) in the electromagnetic calorimeter and is misidentified as an 
electron, and the energies of the remaining jets fluctuate to produce fa. Although the 
probability for a jet to be misidentified as an electron is small, the large cross section for 
QCD multijet events makes this background significant. 

This background is estimated using samples of "good" and "bad" electrons. A "good" 
electron has the quality requirements described in Sec. [IV A 1| , while a "bad" electron has an 
EM cluster with EM fraction > 0.95, Isolation < 0.15, and either if -matrix \ 2 > 250 or track- 
match significance > 10. We assume that the shape of the fa spectrum of the events with a 
bad electron is identical to the fa spectrum of the QCD multijet background. Furthermore, 
with the assumption that the contribution of signal events at low fa is negligible, the bad- 
electron sample can be normalized to the good-electron data in the low- fa region and the 
fa distribution of the bad-electron events can then be extrapolated to the signal region of 
the good-electron sample. 

To estimate the multijet background, we use triggers that do not require fa. Several 
L2 triggers in Run lb meet this requirement, in particular the triggers EM1_ELE_M0N and 
ELE_1_M0N described in Sec. |V|. To avoid biases, we add a condition that the EM object in 
these triggers pass the same L2 requirements as the signal. We then extract two samples 
from these data, based on the electron quality. The fa distribution for the bad-electron 
sample is then normalized to agree with the fa distribution for the good-electron sample at 
low fa (fa < 15 GeV). Figure ^ shows these two distributions. The normalization factor 
Np is calculated as the ratio of the number of bad-electron events to the number of good- 
electron events with < fa < 15 GeV. After imposing the jet selection requirements on 
the events, we find N F = 1.870 ± 0.060 (stat) ± 0.003 (sys). The systematic uncertainty on 
the normalization factor is obtained by varying the range of fa used for the normalization 
procedure from 0-12 GeV to 0-18 GeV. 
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FIG. 6. fa distributions for the good-electron (histogram) and bad-electron (solid circles) 
samples selected from data taken with the EM1_ELE_M0N and ELE_1_M0N triggers (see text). The 
bad-electron sample is normalized to the good electron sample for fa < 15 GeV. 



In the next step, we select two samples from the data taken with the trigger for signal 
events, one containing background and signal ("good" electrons obtained through our se- 
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lection procedure) and the other containing only background events ("bad" electrons). The 
normalization factor Np is then applied to the background sample. Figure [7| shows the 
distributions of fa for the candidates and the estimated QCD multijet background based 
on the bad-electron events after the imposition of jet requirements. 
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FIG. 7. Distributions of fa of the good-electron (sum of signal and background) and 
bad-electron (background only) samples selected from data taken with the trigger used for sig- 
nal events. 



From the above procedure, we estimate 104.3 ± 8.2 (stat) ±9.1 (sys) background events 
for fa > 25 GeV. The systematic uncertainty (8.7%) includes the uncertainty on the nor- 
malization factor (1%), the difference when an alternative method is used to estimate the 
multijet background (5.2%), and the difference for events with fa > 25 GeV when the fa 
region 15-25 GeV is used for normalization (6.9%). In the alternative method, the probabil- 
ity of a jet to be misidentified as an electron is multiplied by the number of multijet events 
that satisfy selection criteria when one of the jets in the event is treated as an electron. 
When more than one jet in an event satisfies the kinematic requirements, all are considered 
in estimating the background from multijet production. 



B. Inherent Background 



The background contribution from processes with similar event topology (i.e., with final- 
state objects identical to those of the signal) is estimated using Monte Carlo events. 



1. W+ > 2 jets 



W+ > 2 jets production is the dominant background to the WW and WZ sig- 
nals. This background is estimated using the Monte Carlo program vecbos [30|, fol- 
lowed by HERWIG for the hadronization of the partons generated in VECBOS and 
then by the detailed simulation of the D0 detector. The cross section from vecbos has 
a large uncertainty, and the generated W+ > 2 jets sample is therefore normalized to 
the candidate event sample after subtraction of the QCD multijet background. To avoid 
the inclusion of WW and WZ events in this normalization procedure, we use only the 
events whose two-jet invariant mass lies outside of the mass peak of the W boson (i.e., 
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Mjj > 50 or Mjj > 110 GeV/c 2 ). Figure |8| shows the two-jet invariant mass distribu- 
tions for data and the estimated background. The normalization factor is found to be 
N v = N VB /(N CSLnd - N QCD /N F ) = 3.41 ± 0.31(stat) ± 0.29(syst), where N VB = 879 cor- 
resonds to the number of vecbos events, N cand = 392 is the number of candidate events 
in the data, and Nqcd = 251 is the number of QCD multijet events outside of the W 
boson mass window. Using this normalization factor, we estimate 279.5 ± 27.2 (stat) ± 
23.8 (sys) W+ > 2 jets events in the candidate sample. The systematic uncertainty is 
due to the normalization of the multijet background (6.9%), uncertainty in the jet energy 
scale (4%), and the difference observed when the range of excluded Mjj is changed to 40- 
120 or 60-100 GeV/c 2 (3%). The cross section multiplied by the branching fraction for 
W+ > 2 jets production, with the W boson decaying to eu, determined with this method 
is 38795/(3.4 x 82.3) = 138.6 ± 14.3 pb (where 38795 is the number of vecbos events gen- 
erated, 3.4 is the normalization factor Ny, and 82.3 pb -1 is the integrated luminosity of 
the data sample), which is consistent with the value (135 pb) given by the VECBOS pro- 
gram. Figure |9] shows distributions in the difference Pt(cv) — pr(jj) and in the separation 
between jets ATZ(jj), which provide sensitive measures for how well background estimates 
describe the jets in the data. The backgrounds from the W+ > 2 jets and QCD multijet 
contributions are seen to agree well with the data. 
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FIG. 8. Dijet invariant mass distribution. The solid circles and the histogram are the candidate 
events and the background events from W+ > 2 jet events and QCD multijet events with a false 
electron, respectively. 



2. tt -> W + W-bb -» eujjX 

Since no limit on the number of jets is applied to retain high efficiency, tt — > W + W~bb — > 
evjjX events contribute to the candidate sample. A sample, simulated using isajet with 
M t = 170 GeV/c 2 , is used to estimate this contribution. We find it to be small, 3.7 ± 0.3 
(stat) ± 1.3 (sys) events. The production cross section for tt events is taken from the D0 
measurement (5.2 ± 1.8 pb) WA - The error in this measurement (35%) is included as a 
systematic uncertainty in our analysis. 
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FIG. 9. (a) Distributions in pT(ev) —pr(jj) before imposition of the mass window on Mjj. (b) 
Distributions for the separation of two jets in r] — (ft space. 



3. WW/WZ -» Tvjj -> evujj 



Since the contribution from WW/WZ — > rz/jj — > euujj is small, and no separate 
simulation of the signal is available, we treat it as background. We use the isajet event 
generator and the detailed detector simulation program to estimate this source. The WW 
and WZ production cross sections are assumed to be 9.5 pb and 2.5 pb, respectively. After 
event selection, we find 0.15 ILqH (stat) ± 0.01 (sys) events. The systematic uncertainty on 
the background estimate is assigned to have the larger value of the asymmetric errors on 
the theoretical cross section (8.4%) p9 |. 



4. ZX -» 



The ZX — > eeX processes can produce events that can be misidentified as signal. These 
events can be included in the candidate sample if one electron goes through a boundary in 
a calorimeter module and is measured as Ifir in the event. From a sample of 10,000 isajet 
ZX — > e + e~X events generated, none survive the selection procedure. The background 
from events of this type is therefore negligible. 
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5. ZX -> t+t-X -» ei/jjX 



The — > rrX processes can also produce events that can be mistaken for signal if, 
due to shower fluctuation, one or two jets from ISR or FSR are observed in the detector. 
From a sample of 10,000 PYTHlA-generated ZX — > rrX events, none survive our selection. 
The background from this source is therefore also negligible. 



IX. RESULTS 



A total of 399 candidate events remain after all selections. The number of events expected 
from SM WW/WZ and from SM background processes are 17.5 ± 3.0 and 387.5 ± 38.1, 



respectively. The transverse mass distribution of the candidate events is shown in Fig. |I0| , 
along with the contributions from background and the SM production of WW/WZ. The 
distributions for data agree well with expectations from background. Table [IV] summarizes 
the number of candidate events, the estimated backgrounds, and SM predictions for the Run 
la and lb data samples. 

TABLE IV. Number of events for backgrounds, data and SM prediction for Run la and Run lb. 





Run la 


Run lb 


Luminosity 


13.7 pb- 1 


82.3 pb^ 1 


Background 


QCD multijet 


12. 2± 2.6 


104.3±12.3 


W+ > 2 jets 


62.2±13.0 


279.5±36.0 


ti — ► evjj + X 


0.87±0.12 


3.7± 1.3 


Total background 


75.5±13.3 


387.5±38.1 


Data 


84 


399 


WW+WZ (SM prediction) 


3.2±0.6 


17.5+3.0 



Figure |TT| shows the pt distributions of the ev system for data, background estimates, 
and SM predictions. We do not observe a statistically significant signal above background. 

Of the 399 events that satisfy the selection criteria, 18 events have pr(eu) > 100 GeV/c. 
The numbers of background and SM events in this pt range are estimated as 18.5 ±1.8 and 
3.2 ± 0.5, respectively. The absence of an excess of events with high pr(ev) excludes large 
deviations from SM couplings. 



A. Limits on Anomalous Couplings Using Minimum pjt 



The WW/WZ production cross section increases, especially at high the coupling 

parameters deviate from the SM values, as shown in Fig. ^|. The p^ distribution for back- 
ground is softer than that of WW/WZ production with anomalous couplings. When events 
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FIG. 10. Transverse-mass distributions of the electron and v (Iftr) system. The solid circles, 
solid histogram, and dotted histogram are, respectively, the candidate events, the background from 
QCD multijet events with false electrons and W+ > 2 jet events, and the expected SM production 
of WW/WZ events scaled up by a factor of ten. 
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FIG. 11. The pt distributions of the ev system from the 1993-1995 (Run lb) data. The solid 
circles are data. The light-shaded histogram is the SM prediction for the background, including 
the dark-shaded histogram, which represents the SM prediction for WW/WZ processes. 



are selected with pf above some large minimum value, almost all background events are 
rejected, but a good fraction of signal with anomalous couplings remains, providing better 
sensitivity to such couplings. This kind of selection eliminates most of SM production, and 
therefore does not have sensitivity to the SM couplings. Moreover the 95% C.L. upper limit 
on the number of signal events (jv 95 « iC - Ij -) can be obtained from the observed number of 
candidate events and the expected background beyond some minimum pjf cutoff, using the 
method described in the report by the Particle Data Group [|]]. To do this, Monte Carlo 
events are generated for pairs of anomalous couplings in grid points of Ak and A. We as- 
sume that the couplings for WWj and WW Z are equal. The expected number of events 
for each pair of anomalous couplings is calculated using the integrated luminosity of the 
data sample, and entered into a two-dimensional density plot with Ak and A as coordinate 
axes. The results are fitted with a two-dimensional parabolic function, and limits on anoma- 
lous couplings are calculated at the 95% C.L. from the intersection of the two-dimensional 
parabolic surface for the predicted number of events with a plane of j\T 95%c,L ' values. The 
resulting contour of constant probability is an ellipse in the Ak — A plane. The numerical 
values for the "one- dimensional" 95% C.L. limits (setting one of the coordinates to zero) are 
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summarized in Table [V| for different minimum values of pjf . 

TABLE V. Limits on A and Ak at the 95% C.L. as a function of minimum pjf , for A = 1.5 
TeV. The number of candidates (A canc j), background (Nbg), and the SM WW/WZ predictions 
(Asm) are also listed. 



(GeV/c) 


-ZVcand 


A BG 


Asm 


A 


Ak 




events 


events 


events 


(Ak = 0) 


(A = 0) 


150 


4 


2.8 


1.9 


(-0.66,0.67) 


(-0.96,1.08) 


160 


1 


2.1 


1.8 


(-0.54,0.54) 


(-0.79,0.89) 


170 





1.5 


0.9 


(-0.52,0.52) 


(-0.76,0.86) 


180 





1.2 


0.2 


(-0.59,0.58) 


(-0.87,0.96) 


190 





0.7 


0.1 


(-0.64,0.64) 


(-0.96,1.05) 


200 





0.3 


0.1 


(-0.74,0.73) 


(-1.13,1.20) 



B. Limits on Anomalous Couplings from the pj. Spectrum 

The limits obtained for some cutoff minimum p^f do not take into account information 
that is available in the full p^ spectrum, and depend on the chosen minimum px value as 
well as on the overall normalization factors for background and predictions for signal. An 
alternative way to proceed is to fit the shape of kinematical distributions that are sensitive 
to anomalous couplings. This usually provides tighter limits, since it uses all the information 
contained in the differential distributions, and it is also less sensitive to overall normalization 
factors. 

As described in Sec. [I], the differential distribution that is most sensitive to anomalous 
couplings is the p% distribution. Our analysis relies on the pt{W — > ev) spectrum rather 
than pt(W — > jj) or pr{Z — > jj) because the resolution on pr(eu) (12.5 GeV/c) is better 
than on Px(jj) (16.7 GeV/c). This is primarily due to the ambiguity in assigning jets to the 
W{Z) boson. 

The differential cross sections have been exploited by previous publications 
P HTT] , P^|HTrj| . p^ ,pD[| for extracting limits on trilinear gauge boson couplings. We use a mod- 
ified fit to the binned p?f distribution to obtain limits, with the modification consisting of 
adding an extra bin in pjf with no observed events, thereby improving the sensitivity to 



anomalous couplings [33 



Based on the number of expected WW/WZ — > evjj events, we choose two 25 GeV/c 
bins between and 50 GeV/c, five 10 GeV/c bins from 50 to 100 GeV/c, two 20 GeV/c bins 
from 100 to 140 GeV/c, one 30 GeV/c bin from 140 to 170 GeV/c, and a single bin from 170 
GeV/c to 500 GeV/c. The cross section for p^ > 500 GeV/c is negligible for any anomalous 
couplings allowed by unitarity. For each bin % of , the probability Pi for observing Ni 
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events is given by the Poisson distribution: 



_ (k + £e^(A, Ak)) n * (h+Cei<7i(XAK)) 
NA 

where £ is the luminosity, and 6j, e^, and and <7j are the expected background, the total 
detection efficiency, and the cross section, respectively, for bin i. Our fast Monte Carlo is 
used to calculate e^a^X, Ak). The joint probability for all bins is the product of the 
individual probabilities Pi, P = n^i" Pi- Since the values C, bi, and q are measured values 
with their respective uncertainties, we assign them Gaussian prior distributions of mean 
fi = 1 and standard deviation o x : 



P' = jGfJfnjGf.dMl 



NA 



where rii = LeiOi is the predicted number of signal events, and Qf n and Qf b are Gaussian 
distributions for the fractions of signal and background events. The integrals are calculated 
using 50 evenly spaced points between ±3 standard deviations. For convenience, the log of 
the likelihood L = logP' is used in the fit and the set of couplings that best describes the 
data is given by the point in the A — Ak plane that maximizes the likelihood given in the 
above equation. 

It is conventional to quote the limits on one coupling when all the others are set to their 
SM values. These "one-dimensional" limits at the 95% C.L., assuming that the WW'y and 



WWZ couplings are equal, are shown in Table [VI]. The limits are more stringent than those 
obtained using the minimum p^f method. 

TABLE VI. Limits on anomalous trilinear gauge boson couplings at the 95% C.L. for three 
values of A obtained using the fit to pjf for data from the Run lb. 



Couplings 




A = 1.0 TeV 


A = 1.5 TeV 


A = 2.0 TeV 


A 7 = X z (Ak 7 = Ak z = 0) 




-0.50,0.53 


-0.42,0.45 


-0.39,0.42 


A/i 7 = Akz (A 7 = Xz = 0) 




-0.66,0.90 


-0.56,0.75 


-0.52,0.70 


A 7 (HISZ) (Ak 7 = 0) 




-0.50,0.53 


-0.42,0.45 


-0.39,0.42 


A/t 7 (HISZ) (A 7 = 0) 




-0.78,1.15 


-0.68,0.98 


-0.63,0.91 


A 7 (SM WWZ) (Ak 7 = 0) 




-1.54,1.58 


-1.53,1.56 




Ak 7 (SM WWZ) (A 7 = 0) 




-2.03,2.45 


-1.79,2.12 




X z (SM WWj) {Ak z = Agf = 


0) 


-0.58,0.62 


-0.49,0.51 


-0.45,0.48 


Ak z (SM WWj) (X z = Agf = 


0) 


-0.86,1.12 


-0.72,0.93 


-0.67,0.87 



We have assumed thus far that the couplings Ak and A for WWZ and WW'y are equal. 
However, this is not the only possibility. Another common assumption leads to the HISZ 
relations p4| . These relations specify Xz, k>z, and gf in terms of the independent variables 



A 7 and k 7 , thereby reducing the number of independent couplings from five to two: Akz 
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~A/« 7 (1 — tan 2 6w), Agf = |A/« 7 / cos 2 9\y, and Xz = A 7 . These one- dimensional limits at 
the 95% CL. are also shown in Table [VI]. 

Since the WW Z and WW'j couplings are independent, it is interesting to find the limits 
on one when the other is set to its SM values. Table [V| includes the one-dimensional limits 
at the 95% CL. for both assumptions: limits on Ak 7 and A 7 when the WWZ couplings 
are assumed to be standard, and limits on Akz and Xz when the WW'j couplings are 
assumed to be standard. These results indicate that our analysis is more sensitive to WWZ 
couplings, as should be expected from the larger overall SM couplings for WWZ than for 
WWj, and that our analysis is complementary to studies of the Wj production process 
which is sensitive only to the WW'j couplings. 



C. Combined Results for Run 1 on WW/WZ — > evjj 



The limits on anomalous couplings presented in this paper are significantly tighter than 



those in our previous publications [11,15]. The primary reason for this is the increase in the 
amount of data (about a factor of six). We can obtain even stronger limits by combining 
the results from Run la and lb. The analysis based on the Run la data is described in 
Refs [lj] , 15| • A summary of the signal and backgrounds for the two analyses [[18 



is given m 



Table IV 



The two analyses can be treated as different experiments. However, because both ex- 
periments used the same detector, there are certain correlated uncertainties, such as the 
uncertainties on the luminosity, lepton reconstruction and identification, and the theoretical 
prediction. Also, the background estimate is common to each experiment. The uncertainty 
on the W/Z — > jj selection efficiency is assumed to be uncorrelated, since we use differ- 
ent cone sizes for jet reconstruction in the two analyses. (This hypothesis does not affect 
the results in any significant way.) The uncertainties for both analyses are summarized in 
Tables [VI I| and [VIII| . Each uncertainty is weighted by the integrated luminosity for the 



respective data sample. Figure [L2| shows the combined plf spectrum. 



TABLE VII. Common systematic uncertainties for Run la and Run lb analyses. 



Source of Uncertainty 

Luminosity 5.4% 

QCD corrections 14% 

Electron and trigger efficiency 1.2% 

Statistics of fast MC 1% 

I$T smearing 5.1% 

Jet energy scale 3.4% 

Total 16% 



To set limits on anomalous couplings, we combine the results of the two analyses by calcu- 
lating a combined likelihood function. The individual uncertainties on signal and background 
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TABLE VIII. Uncorrelated systematic uncertainties for Run la and Run lb analyses. 



Source of Uncertainty Run la Run lb 

ISAJET VS PYTHIA 9% 4% 

Stat, uncertainties of e(W — > jj) 4% 2% 

Parametrization of ej(Tj(A, Ak) 4% 5% 

Total (added in quadrature) 11% 7% 

Background 13% 7% 
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FIG. 12. The pj? spectrum for eujj candidates for the full Run 1 data sample. The solid circles 
are data. The light-shaded histogram is the sum of predictions from the SM and background, and 
the dark-shaded histogram is the SM prediction for WW/WZ processes alone. 



for each analysis are taken into account as in the previous section. Common systematic un- 
certainties are taken into account by introducing a common Gaussian prior distribution for 
the two data samples. 

Combining results from Run la and Run lb yields the 95% C.L. contours of constant 



probability shown in Fig. 13. The one and two dimensional 95% C.L. contour limits (corre- 
sponding to log-likelihood values of 1.92 and 3.00 units below the maximum, respectively) 
are shown as the inner contours, along with the unitarity limits from the ^-matrix, shown 
as the outermost contours. Figure |L3|(a) shows the contour limits when couplings for WWj 
are assumed to be equal to those for WW Z. Figure |13|(b) shows contour limits assuming the 
HISZ relations. In Figs. |l|(c) and |l|(d), SM WWj couplings are assumed and the limits 
are shown for WW Z couplings. Assuming SM WW^ couplings, the U(l) point that corre- 
sponds to the condition in which there is no WW Z couplings (kz = 0, \z = 0, gf = 0) is 
excluded at the 99% C.L. This is direct evidence for the existence of WW Z couplings. These 
limits are slightly stronger than those from the 1993-1995 data alone. The one-dimensional 
95% C.L. limits for four assumptions on the relation between WWj and WW Z couplings: 
(i) Ak = Ak 7 = Akz, A = A 7 = \z, (ii) HISZ relations, (iii) SM WWj couplings, and (iv) 
SM WW Z couplings are listed in Table |T 
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TABLE IX. Limits on anomalous trilinear gauge boson couplings at 95% C.L. from the com- 
bined Run la and Run lb data samples for these values of A. 



Couplings 




A = 1.1) lev 


A = 1.5 lev 


A = z.U lev 


A 7 = \z [Ak^ = Akz = 0) 




—0.42, 0.45 


—0.36, 0.39 


—0.34, 0.36 


Ak^ = Akz (A-v = Az = 0) 




-0.55,0.79 


-0.47,0.63 


-0.43,0.59 


A 7 (HISZ) (Ak 7 = 0) 




-0.42,0.45 


-0.36,0.39 


-0.34,0.36 


Ak 7 (HISZ) (A 7 = 0) 




-0.69, 1.04 


-0.56,0.85 


-0.53,0.78 


A 7 (SM WWZ) (Ak 7 = 0) 




-1.28,1.33 


-1.21,1.25 




Ak 7 (SM WWZ) (A 7 = 0) 




-1.60,2.03 


-1.38,1.70 




A z (SM WW7) {Ak z = Agf = 


0) 


-0.47,0.51 


-0.40,0.43 


-0.37,0.40 


Ak z (SM WWj) (\ z = Agf = 


0) 


-0.74,0.99 


-0.60,0.79 


-0.54,0.72 


Agf (SM WWj) (X z = Ak z = 


0) 


-0.75, 1.06 


-0.64,0.89 


-0.60,0.81 



X. CONCLUSIONS 

We have searched for anomalous WW and WZ production in the eujj decay mode at 
t/s = 1.8 TeV. In a total of 82.3 pb" 1 of data from the 1993-1995 collider run at Fermilab, 
we observe 399 candidate events with an expected background of 387.5±38.1 events. The 
expected number of events from SM WW/WZ production is 17.5 ± 3.0 events for this 
integrated luminosity. The sum of the SM prediction and the background estimates is 
consistent with the observed number of events, indicating that no new physics phenomena is 
seen. Comparing the distributions of the observed events with theoretical predictions, we 
set limits on the WW'j and WWZ anomalous couplings. The limits on anomalous couplings 
are significantly tighter than those using the 1992-1993 data sample. The two results are 
combined to set even tighter limits on the anomalous couplings. With an assumption that 
the WWj and WWZ couplings are equal, we obtain —0.34 < A < 0.36 (with Ak = 0) and 
-0.43 < Ak < 0.59 (with A = 0) at the 95% C.L. for a form factor scale A = 2.0 TeV |35J. 
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FIG. 13. Contour limits on anomalous couplings at the 95% C.L. (two inner curves) and unitary 
constraints (outermost curves), assuming (a) An = Ak 7 = Akz, A = A 7 = Xz; (b) HISZ relations; 
(c) and (d) SM WWj couplings. A = 1.5 TeV is used for all four cases. The U(l) point is the 
expectation with no WW Z couplings. 
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